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Abstract: We consider a time inhomogeneous Cox-Ingersoll-Ross diffusion with positive jumps. We exploit 
a branching property to prove existence of a unique strong solution under a restrictive condition on the jump 
measure. We give Laplace transforms for the transition probabilities, with an interpretation in terms of limits 
of mixtures over Gamma laws. 
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We look at a stochastic differential equation of Cox-Ingersoll-Ross (CIR) type with jumps 

(0) dit = [a{t)- (3{t)it-]dt + f y fi{dt, dy) + a{t) V^t-V dWt , t>0 

J(0,oo) 

which is inhomogeneous in time. Some starting point = ^^o > is fixed. We are interested in 
strong solutions driven by (W, /u) where Brownian motion W and Poisson random measure (PRM) /x 
are assumed independent, such that 

fi{ds,dy) on (0, oo) x (0, oo) has intensity 'd{s)dsv{dy) 

for some fi-finite measure v on (0, oo) with 

(1) / {yM)v{dy) < oo. 







Hence jumps of ^ are positive and summable. The functions a(-),a(-) > correspond to differently 
structured input into equation (o), and there is a time-dependent backdriving force /?(•) > 0. 



The aim of this note is to show that equation (o) has a unique strong solution at least for a restricted 
class of measures ly (theorem 1.1), and to give explicit expressions for the Laplace transforms of the 



transition probabilities of the process (o) (theorem 1.2 and remark 1.3). It wih be seen that input 
in terms of a(-) and input in terms of a(-) induces terms of very simihar structure in the transition 
probabiUties. Thus we extend known results (see Ikeda and Watanabe [IW 89, p. 235]) for the classical 
CIR diffusion 

(2) = [a-PCt]dt + a VCt V dWt , t>0 

(Cox, Ingersoll and Ross [CIR 85] where a, (3, a > are constants) to a time inhomogeneous setting 
with jumps. For the classical CIR diffusion ([2]), it is known that transition probabilities can be 
represented as Poisson mixtures of Gamma laws (combine [IW 89, p. 235] with Barra [B 71, p. 
82], or see Chaleyat-Maurel and Genon-Catalot [CG 06, appendix A3]). The structure of the 
Laplace transform for transition probabilities in (o) according to theorem 1.2 and remark 1.3 below 
indicates additional 'mixtures' (more exactly: weak limits of additional mixtures), in time (due to 
time-dependence of t — > a{t) and t 5(t)) and in space (in the sense of Poisson random measures on 
(0,oo) acting with intensity proportional to v{dy) at times t > 0). 

One motivation for this note is an interest in statistical problems for non time homogeneous processes. 
In a recent work with Kutoyants [HK 09], we have a T-periodic deterministic function of time 
t S{t) in the drift of an SDE, all other coefficients of the SDE being time homogeneous, and need 
as principal assumption - in order to formulate limit theorems - a Harris condition which in order 
to be checked requires explicit knowledge of the transition probabilities corresponding to time steps 
which are integer multiples of T. 

Another motivation for equation (o) is to model the membrane potential in a neuron which - being 
part of an active network of neurons ~ receives from its network both 'continuous' and 'discrete' 
input. Depending on some external stimulus and hence on a time-varying degree of activity of the 
network as a whole, this input is not a function of the membrane potential in the receiving neuron. 
In equation (o), we think of the 'continuous input' part a{s)ds as synaptic input collected through 
the dendrites of the neuron, arriving decayed and delayed at the soma (cf. [BH 06]), whereas the 
'discrete input' part J^^ fj,{dt,dy) models synapses located at the soma itself, or in an immediate 
vicinity (thus near to the recording electrode). There are data sets recording the membrane potential 
in one neuron (belonging to a cortical slice observed in vitro which is stimulated by a potassium bath, 
cf. [H 07, section 3.2] or [J 09, section 5.3]; these data have been recorded by Dr. W. Kilb, Institute 
of Physiology, University of Mainz) which do support a CIR type modelization of the membrane 
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potential. Note that the classical CIR diffusion ([2]) has been used as a model for the membrane 
potential in many papers since Lansky and Lanska [LL 87]. However, only under the effect of quite 
particular experimental conditions (such as electronic stabilisation of the membrane potential in the 
neuron which was the case for the data investigated in [J 09, section 5.3]) real data exhibit a be- 
haviour which can safely be considered as 'stationary'. Inhomogeneity in time (one example in [H 07, 
section 4.6]) seems to be the general situation, and only few work has been done so far in this direction. 



1. Results 

We state our assumptions on the coefficients of SDE (o) in detail. We assume the backdriving force 
/?(•) continuous and nonnegative. In the volatility {t,y) — > o"(i)-\/y+ , the function a{-) is continuous 
and strictly positive on [0,oo). The functions a(-), a(-) modelling time dependence of the input are 
nonnegative, right continuous and locally bounded; they satisfy f^{a + a){t) dt = oo for arbitrarily 
large N. 'fl{ds, dy) denotes the compensator a(s) ds v[dy) of ^{ds, dy), and 



the compensated random measure. Recall that /u and W are independent. Our permanent assumption 
on 1/ is ([T|). Whenever v has infinite total mass (hence ^ an infinite number of small jumps over finite 
time intervals) our construction (strong solutions to (o) through an approximation which relies on a 
'branching property') will require the restrictive condition 



for existence of a strong solution to equation (o) driven by W and fi. Below, at the stages of the 
proofs where ^ comes in, we mention this condition explicitely, working through all other parts of 
the proofs under the permanent assumption ([1]) which ensures summability of small jumps. 

1.1 Theorem: Condition ([3j) implies existence of a unique strong solution for equation (o). 

In equation (o) and in similiar equations with subintervals / of (0, oo) 



IJ,{ds,dy) := n{ds,dy) — Ji{ds,dy) 



(3) 




dit 



ait) - Pit) (t~]dt + I y fiidt, dy) + ait) ^J (t-V dWt , t>0 
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which we shall consider, the process of jumps A,^ = (^t — £,t- I does not depend on ( ^t- I 

V /t>o V J t>o 

Hence, under assumption ([1]) alone, an extension of the Yamada-Watanabe criterion ([YW 71]) on 
the lines of [KS 91, p. 291] guarantees pathwise uniqueness. An extension to more general SDE 
with summable jumps is given in [H 08]. Thanks to pathwise uniqueness, in exact analogy to the 
well known result for classical SDE due to Ikeda and Watanabe ([IW 89, theorem IV. 1.1]), it will 
be sufficient to construct some weak solution to (o), or to similiar equations involving subintervals I 
of (0, oo) as above. We construct weak solutions using a 'branching property' inherent to equation 
(o) (see lemma 2.2). It is this 'branching property' which leads to the restrictive condition ([3]) when 
- for measures having infinite total mass - we pass from subintervals In = [(^njl] with (5„ j to 
/ = (0, 1] and take limits. This will be the topic of section 2; we will finish the proof of theorem 1.1 
in 2.6 below. 



We turn to Laplace transforms of transition probabilities for the process (o). For < s < t < oo and 
A > 0, introduce 

(4) C{s,t) := ^ ^-^exp(^J^ p{u)du^ dv , B{s,t) := exp ^- ^ l3{v)dv^ , 



B{0,t)C{s,t) ' ' ^ ■ B{0,s)C{s,t) ' 



(6) ^.,(A) := jis,t) ^ ^,,^,{X) := [ [1 - e'^-^-'W] K^y) 

Note that ^'s^t(A) is well defined under assumption ([T]) alone. 



1.2 Theorem: Assume ([3]) in order to have a unique strong solution to (o). Assume /?(•) strictly 
positive on [0, oo). Then for < s < t < oo and A > 0, 

is the Laplace transform of the transition probability P {Ct ^ ' \ £,s = u) in the process (o). 



1.3 Remark: In the special case a(-) = a for some constant a > 0, we have 

A ^ 



exp(-/ a{v)^^^t{\)dv\ = f 1 + ) , A>0 
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which is the Laplace transform of the Gamma law T(a,p{s,t)). Hence, in reduction to the special 
case of the classical CIR diffusion ([2]), the Laplace transform given in theorem 1.2 coincides with the 
well known expression given in Ikeda and Watanabe [IW 89, p. 235]. 

See also Overbeck and Ryden [OR 97]. Laplace transforms for transition probabilities in (o) and in 
related equations will be the topic of section 3. The interpretation of transition probabilities in (o) in 
terms of limits of Poisson mixtures of Gamma laws will become clear through lemmata 3.3, 3.7, 3.8, 
3.10. Remark 1.3 will be proved in 3.4+3.5. We will finish the proof of theorem 1.2 in 3.11 below. 



2. Proof of theorem 1.1 

In this section, we prove theorem 1.1. First, a localization argument shows that it is sufficient to prove 
theorem 1.1 in case where in addition to the assumptions stated above 



We shall assume ([7]) throughout this proof section. Second, since pathwise uniqueness holds for 
equation (o) by the Yamada- Watanabe criterion, as in [H 08], we have to prove existence of a 
solution. Our key step is the following lemma. Recall that under the assumptions of theorem 1.1, the 
backdriving force /?(•) > in (o) may vanish. 

2.1 Lemma: For < s < oo and u > deterministic, consider the equation 



(7) 



cj(-) is bounded away from both and oo on [0,oo) 
/?(•), a(-), 5(-) are bounded on [0,oo) . 



(*,s,u) 



Pit) Ct dt + a{t) ^ dWt, t>s,Cs = u 



up to time r := inf{t > s : S,t = 0}, and define ^ = on [[r, oo[[. Write 




for the unique strong solution. Then r is finite a.s., and we have for fixed T G (0,oo) 




) 



where C depends on upper and lower bounds for C7^(-) according to ([7]), and not on s, u, T. 
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Proof: 1) The solution of (*, s, u) is unique and strong in virtue of the Yamada-Watanabe criterion 
(see Yamada and Watanabe [YW 71], Ikeda and Watanabe [IW 89, theorem IV. 1.1], Karatzas and 
Shreve [KS 91, p. 291]; note that on every [[0, p^]], pi := mf{t : S^t < j}, we may use standard results 
under Lipschitz assumptions). We fix s = 0; the Markov property implies that is sufficient to prove 
the assertions of lemma 2.1 in this case. 

2) It is sufficient to consider case a{-) = 1 only: the deterministic time change 

A{t) := [ a^{v)dv , r(t) := mf{v > : A{v) > t} 
Jo 

associates to ^ := {S,T(t))t>o a Brownian motion W such that (*,s,u) is transformed to 

(8) dl = -mldt + ^ dWt, t>0, Co = « 

where /3(t) = f3{T{t)) is continuous, nonnegative and bounded; this follows from 

rr{t) rrit) I — 

Cr(t)-?r(5) = - / P{v)i^dv + / a{v)^JitdW^. 

Jt{s) •^I'is) 

By assumption ^ on upper and lower bounds for (t(-), we have positive constants c, c such that 
c{t — s) < T[t) — t{s) < c{t — s). The time change t — > T{t) being deterministic, it is sufficient to prove 
the assertion of lemma 2.1 in the form dH]). Hence we assume a{-) = 1 from now on. 

3) Consider first the unique strong solution fh for the problem without drift 

(9) dfht = \J fhf dWt , t > , rho = u 

and write p for the time where fh enters 0. We have p < oo P-almost surely by [IW 89, p. 237]. fh 
is a local martingale before time p, and a nonnegative supermartingale on [0, oo). Introduce stopping 
times p^^^ < p and stopped processes fh^^^ 

:= inf{t >0:fht^ [-^,£]} , := (m,^^w)t>o 

for i large enough. Then fh^^^ is a uniformly integrable martingale, and ([9]) yields 



(10) E ((r^W - fhi;^)r) = E ((r^W - mj^^)^,,.,) 
for all T. By p^^^ | p as ^ — > oo and by monotone convergence 



E I / mf ^ dt] < u-T 







Wn2 



E sup (mt) = sup E sup (m). ) 
\0<t<TAp J eeJN \0<t<T J 

Since fh is absorbed in at time p, this shows in combination with pU|) 

(11) e{ sup (fhtf] < sup 2 I u^ + e{ sup (mf^ -m{,^^)M I < C {u^ + u ■ T) 

\o<t<T I £elN \ \o<t<T I I 



for some constant C, using [IW 89, p. 110]. On the right hand side of ()lip . the term which is hnear 
in u arises as a consequence of the square root form of the diffusion coefficient. 

4) For the process ^ solving equation ([8]), define r := inf{t > : = 0}. /?(•) being nonnegative, 
we may apply [KS 91, pp. 291+293] which does not need Lipschitz conditions, and compare pathwise 
the solution of ([5]) to the solution of @: we obtain first 

< fnt for all t G [[0, r[[, P-almost surely 

and second 

Ct < rhf for all t £ [0, oo) 

since S^t is absorbed in at time r whereas fh stays nonnegative. Thus step 3) implies t < oo P-almost 
surely, and we get in combination with ()lip 

E \ sup i^tf ] < E\ snp {mtf ) < C {u^ + u-T) . 

\0<t<T J \0<i<r J 

This gives the bound in the asserted form for equation ([5]). The proof is finished. □ 
Similiar to Protter [P 05, p. 250], we will consider for p = 1, 2 and for < T < oo norms 



(12) mi^T := [e{ sup \Xs\n\ 

for cadlag processes X = {Xt)t>o- Such norms with p = 2 appear in lemma 2.1. 



2.2 Lemma: For s > 0, i = 1,2 and ui,U2 > 0, consider strong solutions 

^is,Ui) q£ equations (*, s, Ui) driven by W^"^^ 
where W^'^\ i = 1,2, are independent: then the sum 

is a weak solution ^ to equation (*, s, ui + U2)- 



Proof: This is the well-known 'branching property' of the CIR diffusion. Write 
prepare a third standard Brownian motion l^(^) independent of W^^\ i = 1,2. Then 



^is,u,)^ and 



dWt 



E 



(0) 



is again a standard Brownian motion. Equations {*,s,Ui) driven by W^^\ i = 1,2, show that 
^ = S,^^^ + is a solution for equation (*, s, ui + U2) driven by W. □ 



2.3 Lemma: For < a < 6 < 00 arbitrary, fixed, consider the equation 

(+, (a, b]) dit = -m it- dt + [ y fi{dt, dy) + a{t) J (1-)+ dWt , t > , ^0 = 

J(aM 



driven by W and /i, both independent. 



a) Pathwise uniqueness holds for the SDE (+, (a, 6]). 

b) Arrange the point masses of /i on (0, 00) x (a, b] in increasing order of time 

00 

IJ.{ - n (0, 00) X (a, b]) =: ^{Ti,Yi) '^^lere < Ti < T2 < . . . and lim Ti 

1=1 

prepare Brownian motions W^^\ i > 1, independent and independent of and strong solutions 

^{T^,Y^) for SDE {+,Ti,Yi) driven by i > 1 . 



Then 



is a weak solution to (+, (a, 6]). 



^ = ^«'''] := ^^m.^i) 1 



[[r„oo[| 



1=1 



c) Equation (+, (a, 6]) has a unique strong solution ^ = F(Wi /i). 

d) Using the norms defined in (fT2]) . the solution of (+, (a, 6]) satisfies 



(13) 



*1T 



J{a,b] 



(14) 



*2T 



< CT{1+Tf u 



J{a,b\ 



for suitable constants C which only depend on upper and lower bounds for cr(-) and on upper bounds 
for a(-), according to ([7]). 



Proof: 1) Extend the nonnegative processes of lemma 2.1 with starting time s > to the full 
time interval [0, 00) by Sjf'^^ = for < t < s. For any collection of points {si,Ui) where Ui > and 
< si < S2 < . . . such that lim Sj = cxd, write ^(**'"») = F{W^, Si,Ui) for strong solutions of (+, Si,Ui) 



driven by independent Brownian motions W^. Fix < T < oo. As a consequence of lemma 2.1, and 
by the norm properties of ()12|) . we have 



(15) 



i=l 



< 



*1T 



i=l 



< ^C{l + T) Y^Ju, + u] 



*2T 



1=1 



2) With probabiUty 1 we can arrange the points of ^{dt,dy) on (0, cxd) x (a, 6] as 

(Ti,Yi) , i > 1 , such that < Ti < T2 < . . . and hm Tj = 00 . 

j— >oo 

In this representation, (Yi)i are iid with values in (a, 6] and with law The sequences 

(Tj)j and are independent. Similiar representations have been used in [HJ 93]. The number 

Nt ■= ^{ (0, T] X (a, b] ) of points (Tj, 1^) with Tj < T is Poisson distributed with parameter 

v{{a,b]) [ a{s)ds < C T i^{{a,b]) 
Jo 

by definition of Poisson random measure with intensity a{s)dsv{dy), and by assumption on a(-). 

3) For PRM ^ of step 2) independent of the BM's VF*, i > 1 of step 1), we may read inequality p5|) 
with p = 2 and both sides squared as an inequality for conditional expectations given fi: 



(16) 

where 



E 



Nt 

sup E^P'"' 



< Cil + T) 



Nt 



.1=1 



= f^{/u((0,t]®S) : t > 0,5 G ^(iR)} 



denotes the cr-field generated by /i and where 

Nt ^ 



~{ ^ Jo J{a,b] 



y2 fi{ds,dy) . 



Using the elementary inequality (a + b) < 2(a + b ), we get 
2 

< 2 



.i=l 



T 



J(a,b] 



Vy + y'^ fJ-ids,dy) 



- 2 




+ 


/7 

Jo J{c 



Vy + y^ Hds, dy) 



Here the second term on the right hand side is deterministic since 'fl{ds,dy) = a{s)dsv{dy); using 

bounds for a(-) according to ([7|) and Jensen inequality, this term is smaller than 

2 



CT^ {v{{a,h])y 



/ Vy + 

J{aM 



u{dy) 



< CT^iyi{a,b]) / {y + y^),.{dy) . 

J(a,b] 



l(a,b] '^{ia,b]) 

The first term on the right hand side is the square of a martingale and thus has expectation 



{y + y^)Jl{ds,dy) < CTu{{a,b]) / {y + y^)u{dy) 

J(a,b] J(a,b] 
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where we use again ([7]). Putting all this together, we obtain 



(17) 



E 



Nt 

sup EeP'"'^ 



< CT{l+Tfiy{{aM) [ {y + y'')v{dy) 

Jia,b] 



for some constant C which only depends on upper and lower bounds for a{-) (through (jl6p and lemma 
2.1) and upper bounds for a(-). The right hand side of ()17p corresponds to ()14|] . 

4) To prove (fT3l) , we read (fTSl) with p = 1 conditionally on /i as 



(18) 



E supJ^^f'^'M-^J < ^/c(^T^) ^ . 



Again, by the bounds on a(-), we have 

e( r [ Vy + y^ Kds,dy)] = r [ Vy + y^Kds,dy) < ct [ + y^ ^{dy) 

\JG J(a,b] J Jo J{a,b] J {a,b] 



and thus 
(19) 

where the right hand side of corresponds to ([13]). 



£;f sup Vef"'''^! < CTVT+T [ 7^T^z.(dy) 



5) We show that with notations as above, the process 

provides a weak solution to equation (+, (a, 6]). To see this, prepare as in the proof of lemma 2.2 an 
additional standard Brownian motion W^^^ independent of W^^\ i G IN, and consider 



dWt 



1, 



E 



dH^/') + 1, 



Then W is again a standard Brownian motion, and we have 



(0) 



d^t = -Pit) Ct- dt + y fi{dt, dy) + a{t) {^t- )+ dWt , t > , = 

J{a,b] 

by definition of the processes 

^{T„Y^)^ i £ Note that in restriction to every interval [[0,T„,[[, the 

sums reduce to finite sums. 

6) We finish the proof of lemma 2.1. First, pathwise uniqueness holds for solutions of equation 
(+, (a, b]) by Yamada-Watanabe, as in [H 08]. Second, we have constructed a weak solution to (+, (a, b]) 
in the preceding step 3). Hence, as in Ikeda and Watanabe ([IW 89, theorem IV. 1.1]), there is a unique 
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strong solution for equation (+, (a, b]) which has the form of a functional F{W, ^) of the driving pair 
{W,fj.). We have to justify the last assertion: 

Inspecting the proof [IW 89, p. 163-166] for classical SDE that pathwise uniqueness and existence 
of a weak solution together imply existence of a strong solution, functional of the driving Brownian 
path ([IW 89, theorem IV. 1.1]), we see that this proof carries over to SDE with jumps where the 
solution process is cadlag. In fact, the proof is based only on the following: i) cadlag driving 
processes with stationary and independent increments; ii) transition probabilites between path 
spaces which are Polish. So we can replace C used there (as path space for the solutions) by 
D, and Cq = {a £ C : q(0) = 0} used there (as path space for the driving process) by CqxDq 
since Dq = {a £ D : a{0) = 0} is a path space for ^ J^^ y fj.{ds, dy)^ under assumption ([T]). □ 

Now we turn to the problem 



For measures v where i^{{0, 1]) is finite, equation (+, (0, 1]) has finitely many jumps in finite time 
intervals, thus we can proceed as in lemma 2.3 to construct a weak solution and - thanks to path- 
wise uniqueness - a unique strong solution ^ = F(W,fj,). It remains to consider measures v having 
z^((0, 1]) = oo. If we think of a tending to in lemma 2.3 and compare the structure of bounds p3p 
and (|14|) for a | 0, the restrictive condition ^ 






(0,oo) 



seems unavoidable in this case. We start with a weak solution for (-I-, (0, 1]) under ([3j). 



2.4 Lemma: Assume ([3]) and i^{{0, 1]) = oo. Fix some sequence 5n I such that 



(20) 



< 




—n 



n > 1 



and 6o = 1- Prepare BM's W^^\ n > 1, independent and independent from /i, together with 



^S„,s„-i] strong solution to (-|-, {6n, 



'n-l 



]) driven by (VF("),/x) 



for n > 1. Define for > 1 



N 




n=l 
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Then the paths of f^^) converge uniformly on compacts as A'" — > cx) to the paths of ^ = 



(21) C^O'i] := supp) = 



n=l 



where the process defined by ()21|) is cadlag, provides a weak solution to the problem (+, (0, 1]) , and 
satisfies the bound 

(22) ^O'l] < ctVTTt [ Vy + y^ ^{dy) 

i(0,l] 

for some constant C which depends only on bounds for the functions in ([7]). 



Proof: We assume that the total mass z^((0, 1]) equals +oo, and fix < T < cxd. By choice (pOj) of 
{5n)n-, by the bound (fT3l) in lemma 2.3 combined with Chebychev inequality, we have 

P\ sup g'^"'-'"-!] > 2~" ) < CT^TTT 2"'^ 
\o<t<T y 

for all n > 1. Using Borel Cantelli, we see that P-almost surely the paths of 

TV 



n=l 



converge uniformly on [0, T] as ^ oo to a finite limit 

i := supp) . 
N>1 

A uniform limit of cadlag processes is again cadlag. By monotone convergence, and applying the 
bounds ()13p of lemma 2.3 to all processes N > 1, we obtain for the limit ()2ip the bound 



E 



I sup l] < ctVi + t [ Vy + y^ <dy) 

\0<t<T J 7(0,1] 



thanks to assumption This is the bound 

2) We prove that the limit in ()2ip is indeed a weak solution to equation (+, (0, 1]). Due to the bound 
([22]) . the process 



rt I 

^ 1 



is a well defined martingale whose angle brackett 



°o ft ft 

'0 Jo 



n=l 
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is integrable at every fixed time t < oo. Prepare as in step 4) of tlie proof of lemma 2.3 another 
standard Brownian motion W^^^ independent of W^"'\ n G IN, and of /u, and consider 



■= %->0}E 



2~((5„,5„_i] 



n=l \ 

Then W is again a standard Brownian motion, and we do have 

d6 = -Pit) I- dt + [ y fi{dt, dy) + a{t) ^(6-)+ c^VP^t , t > , = 

J(o,i] 

by definition of the processes ^i^^^^^ri-i] ^ ^ here we exploit convergence ^^^-^ | ^ a.s. uniformly on 
compacts as ^ oo, proved in step 1), to show convergence of drift terms 

Pit)^^^^dt — y [ Pit)^t-dt , iV^oo 
'0 Jo 

and angle bracketts of martingale terms 

/^l <^{t) dWt - Y.J^<y{t)^^y-^^ dW^-^^ , iV^oo 

a.s. uniformly on compact time intervals. Since /i induces positive jumps only which are summable, 
it is obvious that 



yfi{ds,dy) — ' yKds,dy) , N ^ oo 

/O ^(5]v,l] -'0 J(o,i] 

a.s. uniformly on compact time intervals. The three last assertions together prove that the limit 

process ^ in (f2T]) provides indeed a solution to the problem (+, (0, 1]) . □ 



2.5 Lemma: Under assumption ([3]), for measures v with i^((0, 1]) = oo, 

a) pathwise unicity holds for equation (+, (0, 1]); 

b) equation (+, (0, 1]) has a unique strong solution ^ = F{W,^). 

Proof: Exactly as in step 6) in the proof of lemma 2.3, based on existence of a weak solution which 
has been established in lemma 2.4. □ 

2.6 Proof of theorem 1.1: 1) In order to solve equation (o) with starting point xq > 0, introduce 
BM's W^^\ i = 1,2,3, independent and independent of /i, together with the following processes 
^» =F(TyW,/x): 
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i) ^^^^ = F(W^^\ fi) is a strong solution to equation (+,(0,1]) according to lemma 2.5, driven by 
W := W^^^ and the small jumps of /x; 

ii) ^(^^ = F{W^'^\ fx) is a strong solution to equation 

(23) dCt = -m it- dt + [ y fiidt, dy) + (j(i) ^^^^t , i > , Co = 

J(l,oo) 

driven by := W^'^^ and the big jumps of (since by assumption there are at most finitely many big 
jumps over finite time intervals, this strong solution is obtained as in steps 2)+5)+6) of lemma 2.3; 
an analogue of (jl6p holds, whereas no finite moments and hence no moment bounds corresponding to 
(fT3l) or (fH|) are available for ^^^^ = F{W^'^\^) under our assumptions on v); 

111 is a strong solution to the 'continuous' problem 



(24) d^t = [a{t) - Pit) it] dt + ait) V(^ dWt , t > , io = xo > 

with xq > the starting point of equation (o) (in case xq = 0, the process (I24p is identically 0; in case 
xq > 0, a unique strong solution to (pl|) exists as in lemma 2.1). 

2) Prom the processes defined in step 1), we construct a weak solution ii,W) for equation (o) with 
starting point xq by 



i=i \ st- 



i=l 



Since pathwise uniqueness holds for equation (o) (extension of Yamada-Watanabe as in [H 08]) we 
conclude as in step 6) of the proof of lemma 2.3 (extension of [IW 89, theorem IV. 1.1] to SDE with 
jumps) that there is a unique strong solution i = FiW, fi, xq) for equation (o), functional of the 
driving pair iW, fi) and of the starting point xq. Theorem 1.1 is now proved. □ 



3. Proof of theorem 1.2 

We give the proof of theorem 1.2 through a series of lemmata. Again, a localization argument shows 
that it is sufficient to consider the case where the assumptions of 1.2 are strengthened to 

/?(•) and o"(-) are bounded away from both and oo on [0, cxd) 
(25) ^ 

a(-), «(•) are bounded on [0,cxd) . 
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We shall assume (1250 throughout this proof section. We start with asymptotics of (jlJ+dS]) under 
condition ([25]), and use the notation at x ht in the sense of < liminf ^ < limsup ^ < cxd. 



3.1 Lemma: Under ([25D, for < s < t < cx) and 7(s,t), p{s,t), B{s,t) defined in ([3])+(l5]): 

2 1 

(26) 7(s, t) ~ as t i s , 

2 1 

(27) p{s,t) ~ ~Yr^7Z~ as t I s , 

(28) p{s,t) X 1 as t T oo , 

(29) 7(s,t) X -B(s,t) as t T oo . 

Under (I25p . for every choice of < s < oo, B{s,t) decreases exponentially fast as t — > oo. 



Proof: For < s <t and t decreasing to s, we have for C{s,t) and B{s,t) defined in (jlj) 



t - s 



C{s,t) 



B{0,s) 



2 exp y P{v)dv 



B{0,s) 



as t i s , 



which gives first ([26]) by definition of 7(s,t) in ([5]), and second (f27|) for p{s,t) 



B{s,t) ■ 



For s fixed and t increasing to oo, assumption ([25]) guarantees 



B(0,t)C(s,t) 



exp 



P{r)dr I 



as t ^ oo 



which gives first ([28]) by definition of p(s,t) in ([5]), and second ([29]) for 7(5, t) = B {s , t)p{s , t) . □ 



3.2 Lemma: Consider < s < t < 00. For A ^s,ti^) defined in ([6]) we have 
and for all A > 

(31) lim^',,t(A) = A, lim^,,t(A) = 0. 

tis tfoo 

Moreover, the family (^s,i)o<s<t<oo ^ semigroup property under functional iteration: 

(32) o *t2,t3 = *ti,i3 on [0, 00), for < ti < f2 < t3 < 00 . 



Proof: Calculating derivatives with respect to A in 
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at A = O"*" gives (I30p . The limits (I3ip are an immediate consequence of lemma 3.1. We turn to the 
iteration property. Fix < ti < ^2 < ^3 < oo. First, 

(33) ^'''■'-'^^fe-'^>- mJl5':)cfe.*3) 

follows from definition of p(ti,t2) and 7(^2,^3) since 

1^1 _ 1 C(ti,t2)+C(t2,i3) 



B{0,t2)C{ti,t2) B{{)M)C{t2,h) Bi0,t2) C{ti,t2)C{t2,t3) 

where the numerator is by definition additive. Using (133]) . the definition oi p{., .) and 7(., .) gives 

(34) Piti,t2)pit2,h) ^ 7ih,hhit2,h) ^ 

^ ' P{tl,t2)+l{t2,h) ^^1'^^' p{h,t2)+-f{t2,t3) '^^''^ 



Using (|33|) + (jM|) . we calculate 



P{tl.t2)+l{t2M) pit,,t,)+X 

iitut2Mt2,h) A 



Pitl,t2)p{t2,h) + \{p{tut2) + 7(^2, is)) 



which is the iteration property (|32|) . In this last part of the proof, we follow closely the representation 
given in Hammer [H 06] for the classical CIR diffusion. Lemma 3.2 is proved. □ 



The next step is to identify a semigroup {Hs^t{-, '))o<s<t<oo transition probabilities on the state 
space ([0, cxd), ;B([0, 00)) of (o) with the property 

(35) Hs,t{u, B) = P (^^l"'"'^ e , < s <t < 00, u>0, B e B{[0,oo)) 

where ^^^'^^ denotes the solution to equation {*,s,u) in lemma 2.1. For the classical CIR diffusion, 
Poisson mixture representations of the transition probabilities have been considered in Chaleyat- 
Maurel and Genon-Catalot [CG 06, (78)+(82)]; see also Going-Jaeschke and Yor [GY 03] for some 
context. 

3.3 Lemma : a) For < s < t < 00, introduce a transition probability Hs^t{', ■) 

^^s,t(y,-) := e-^^(^'*)eo(-)+E^ r(m,p(.,t)) , y>0 

m=l 
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on ([0, oo), S([0, oo)). For y > 0, Hs^tivr) is a Poisson mixture of Gamma laws; for y = 0, -ffs,t(0, ■ 
reduces to Dirac measure eo sitting in 0. For all y > 0, the Laplace transform 



/•oo 

/ e-^"" Hs,t{y,dx) ,A>0 
Jo 



of Hg^tiy, •) is given by Barra [B 71, Ch. VII. 1] as 

¥'s,t(y;A) = exp ^-y 7(s,f) -^-^-p^^ = exp (-y ^'^,t(A)) . 

b) For fixed s and y > 0, under ([25]) . mean and variance of Hs^t{y-, •) 

(36) / xHs,tiy,dx) = yB{s,t), / {x - yB{s,t)f Hs,t{y,dx) = y- 

Jo Jo 

vanish exponentially fast when t tends to oo. 



p{s,t) 



c) {Hg^ti', •))o<s<i<oo is a semigroup of transition probabilities on ([0, oo), ;B([0, oo))). 

d) Write for the class of C°°-functions having compact support contained in (0, oo). In restriction 
to ((0, oo),,B((0, oo))), (-f^s,t(') •))o<s<t<oo acts as a submarkovian semigroup with the property 

gy/-/ ^^^ = (A/)(y) := -myfiy) + ^yf"iy) , / e C y > o 

tls t — S Z 

for every s. 

e) As a consequence of d), the semigroup {Hg^ti', •))s<t on ([0, oo), i3([0, oo))) corresponds to the 
system of SDE {(*, s, y) : < s < oo, y > 0} with absorption at considered in lemma 2.1. 

Proof: 1) Assertion a) follows from [B 71, p. 82]; e) is immediate from d); b) comes from derivatives 
of (ps,tiy', •) at 0+, using lemma 3.1 and ([50]) in lemma 3.2. We have to prove c)+d). 

2) To A > we associate z = e~^, < z < I, and define for < s < t < oo and y > 

9s,t{y;z) := ips,tiy;-'^og{z)) = / z"" Hs,t{y,dx) , o<z<i. 

J[0,oo) 

Then the function gs,t{y] ■) is strictly increasing on (0, 1], takes the value gs,t{y] 1) = 1 at z = 1, and 
satisfies by part a) of the lemma 

5s,t(y; •) = (5s,t(l; for O < s < t < oo and y > . 

It is easy to see that the set of functions 

y,,t(l;-) : (0,1] ^(0,1] , < s < t < oo 
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has a functional iteration property 

(37) 541,^2(1; •) o 5*2,43(1; •)= S'ti,t3(l;-) for < ti < t2 < ts < 00 

similiar to the iterates of the generating function in classical branching process theory (cf. Athreya 
and Ney [AN 72]): this is a consequence of ()32p in lemma 3.2 since by definition 

From (I3ip in lemma 3.2 we have immediately 

(38) Ymi gs,t{'^;z) = z , lim 5s,t(l; 2:) = 1 , 0<z<l. 

3) We prove the semigroup property of {Hs^t{-, ■))o<s<t<oo d^' °°)' '^([0' thus c): consider 
arbitrary points < < t2 < ^3 < 00, then for z £ (0, 1] 

/ / Ht^,t2iy,du)Ht2^t3iu,dx)z^ = / Ht^^t2iy,du)gt2^t3{u;z) 

J [0,00) J [0,00) J[0,oo) 

-f^ti,t2(y,dn) [5(t2,t3(i;^)]" = 9h,t2{y;9t2,t3{i;z)) 

[0,00) 

= [9h,t2{'^;gt2,t3{i;z))]y = [5^^,43(1; z)]^' = gt,,t3{y;z) = / Ht,fy{y,dx)z'' . 

J [0,00) 

4) To determine the generator of {Hg^ti-, ■))s<t at t = s when the semigroup is restricted (0,oo), we 
give a heuristic argument. Consider y > 0. By (I36|) combined with (I27p and by definition of B{s,t) 
in dl]), the laws Hs^t{y-, •) are concentrated as t | s on small neighbourhoods around their mean value 
y B{s,t), which tends to y. Thus we have for / £ as t | s 



iH,,tf - f) (y) = J Hs,tiy, dx) [f{x) - f{y)] 

= [f{yB{s,t))-f{y)] + f'{yB{s,t)) I Hs,t{y,dx)[x-yB{s,t)] 
+ ^f"{yBis,t)) J Hs,t{y,dx)[x-yB{s,t)f + ... . 



On the right hand side of the last equation, the second term is by (l36]) . the third one equals 
y^^f"{yB{s,t)) ~ {t-s)^^yf"{y) as U ^ 

by ([36l) + (f27|l . and the first term behaves as 

f'{y)y[B{s,t)-l\ {t-s)(3{s)yf{y) as t[s. 

Together, this shows for < s < t < cxd and for y > 

lim ^f~U y) = (A/)(y) 

tls t — S 
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with As as defined in d). This gives d), and concludes the proof of lemma 3.3. 



□ 



3.4 Lemma: For < s < t < oo, we have 

This is the Laplace transform of an exponential law with parameter p{s,t). 



Proof: Taking first the log of both sides and deriving then with respect to s for fixed t, we obtain 

2 



2 / \ 

on the left hand side and on the right hand side the same derivative —4-^ ^s^t(A) ; when calculating 



the derivative on the right hand side, we make use of 



3.5 Proof of remark 1.3: This is an immediate consequence of lemma 3.4. □ 

The next step is to specify the semigroup (Js,t(-, ■))o<s<t<oo ^ diffusion ([2^ 

d^t = Ht) - (3{t) 6] dt + a{t) 7(6)+ dWt, t > 
with time dependent coefficients. We start with 

3.6 Lemma: For given a(-) nonnegative, right continuous and bounded, there is a family of laws 
iIs,t)Q<:s^t^ao ([0) oo), i3([0, oo))) with the 'skew convolution' (Li [L 93]) property 

(39) It^^t2Ht2,tz * It2,t'j, = Itify for all < ti < t2 < is < oo 
relative to the semigroup (-f^s,t)o<s<j<oo lemma 3.3, and with Laplace transforms 

(40) ^ e-^"" Is,t{dx) = exp ^- ^ ^-^aiv) ^'t,,t(A) dv^ , 0<s<t<oo,A>0. 

Proof: 1) Note first that for constant a > and ti < t2 < t^, the law 

J[0,oo) 
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has by 3.3 and 1.3 and (I32p the Laplace transform 

2) Fix < s < t < cxD. For the grid of points tnj '■= j S ]Nq, introduce step functions 

oo 

i=i 

which - by right continuity of a(-) - converge to a(-) as n oo. For j > 1 such that intervals 
[tn,j-i,tn,j) and {s,t] intersect, write Qnj for the law on [0,oo) having Laplace transform 



A — ^ exp - / —^aniv)^y,tWdv 



as in step 1). Convolution of the Qn,j defines a law on [0, cxd) whose Laplace transform is 



(41) A ^ exp [- 1^ ^ an{v) ^,,t(A) dv 

Now the right hand side of (j40p viewed as function of A approaches the limit 1 as A | 0, and is the 
pointwise limit of Laplace transforms ()4ip as n ^ cxd. Hence the function ()40p itself (Feller [F 71, p. 
431]) is Laplace transform of some probability measure Is^t on [0, cxd). 

3) Rephrase steps l)+2) as follows: Is,t arises as weak limit of convolutions 

r ( a{rn,i) , p{rnfl, rn,i) ) -f^r„,i,i * . . . * r ( a{rnj) , p(r„j_i, r„j) ) i^r„j,i * • • • 

along grids s = r„^o < '^n,i < • • • < '^ra,i(n)-i < fn,l(n) = i in (s,t] with mesh tending to 0, as n ^ oo, 
or equivalently as weak limit of laws 



n — > oo 



(42) £ |^^^K.„c/„.,) 

where at every stage n of the asymptotics, Un,j ~ F ( a{rnj) , p{rnj^i, rnj) ) are independent random 
variables and ^(^ri,j,u„j) independent processes solving {*,rn,j,Unj) of lemma 2.1, for 1 < j < i{n). 
Note also that with high probability, the random variable Unj takes its values in small neighbourhoods 
of 2 ' Ci{i"n,j) {fnj — ^nj-i), as a conscquence of (P7|) in lemma 3.1. From interpretation ([1^ of 
Is^t, tbe skew convolution property ([39]) is obvious. □ 



The interpretation of Ig^t in terms of limits of laws (j42p reappears in the following lemma: 
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3.7 Lemma: For given a(-) nonnegative, right continuous and bounded, the semigroup of transition 
probabiUties associated to a CIR diffusion with time dependent coefficients 



(43) d^t 



-ait)- Pit) 



dt + ait) dWt , t>0 



2 

(note the particular form of the input term in comparison to (j24|) ) is given by 

iJs,ti-r))o<s<t<oo ' Js,tiyr) = Hs,tiy,-)*Is,t 
where (/s,t)^^^ is associated to a(-) by lemma 3.6, and iHs^t)^^^ is the semigroup of lemma 3.3. 

Proof: For < s < t < oo and y > 0, define Js^tiv-, ) '■= Hs^tiu,-) * h,t- Combining ([5^ with 
the branching property according to lemma 2.2, we check that (Js,i(-, ■))o<s<t<oo ^ semigroup of 
transition probabilities on ([0, oo), ;B([0, oo))). We have to verify that this semigroup corresponds to 
the process (|43|) . As in lemma 3.3, we give a heuristic argument. Since Js,ti', ■) is a convolution, 



iJs,tf - f) iy) = J Js,tiy, dx) [fix) - fiy)\ 

Hs,tiy, dv) I [fiv) - fiy)] + [ IsAdu) [fiv + u) - fiv)] 



Consider the first expression on the right hand side: in case y = 0, this term is (since Hg^tiO, ■) is the 
Dirac measure at 0); in case y > 0, it behaves for t | s as (t — s)iAsf)iy) in the notation of lemma 
3.3 d). We exploit the moment structure of I^,*: from ([l0]) + ([30]) + ([27|) we have 

Eils,t) = ^-Y^aiv)Biv,t)dv ~ {t-s)^-^ais) as 

whereas Var ils,t) is proportional to (t — s)^. In case y = 0, the second expression above is 

Is^dn) [fiy + u)- fiy)] ~ it-s)^ ais) fiy) + O ((t - s)') , Us. 

In case y > 0, the second expression above is 

IsAdn) [fiv + u)- fiv)] ~ it-s)^ ais) fiv) + O ((t - s)^) , Us 

where values v selected by Hg^y^dv) are concentrated on small neighbourhoods of their expectation 
y Bis,t), cf. (f36]) . hence close to y as t | s. Hence in both cases y = (where by lemma 3.3 (^s/)(0) 
equals 0) and y > we end up with 

= {(A/)(y) + ^a(s)/'(y)} , y>0 
which corresponds to equation (j43p . □ 
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The next step is to specify the semigroup ( •) ) associated with equation 

V ' /0<s<t<oo 

(44) d6 = -Pit) Ct- dt + [ y tiidt, dy) + a{t) y^{Ct-)+ dWt , t>0 

J(0,oo) 

where only the 'discrete' part of the input into equation (o) appears. We start with 



3.8 Lemma: Let satisfy the restrictive condition Fix < ti < t2 < t^, and write fi{ds, dy) for 
the restriction of the random measure fi{ds,dy) to (ti,t2]x{0,oo): then 

fl{ds,dy) is PRM on (ti, x (0, oo) with intensity i>{ds,dy) ■.= 'd{s)l^i^^i^^{s)ds v{dy) . 

Represent /x in form ^{Ti.Yi) ^-^^d associate to the point masses of Ji sitting at (Tj,li) processes 

^^i"^"^'^) which are strong solutions to (*,Tj,li) according to lemma 2.1, driven by independent 
BM's and hence independent, conditionally on /Z. Then thanks to ([3|), the sum 

oo „ 



(45) E4 =■ / i\rms,dy) 



i=l 

is finite a.s. and has the Laplace transform 



(ti,t2]x(0,oo) 



(46) A — > exp<^ - / i-g-s'^-.'sW v{ds,dy) 

y i(ti,t2]x(o,oo) L -I 

or with notation of ([6]) 

'■t2 



A — > expj-^ a(s) ^s,f3(A) (is I , A>0. 



Proof: 1) Consider first some a > and some finite collection of points (sj,Uj) in (ti, x [oj oo), 
1 < i < £. Associate to the points {si,Ui) processes fej'""'^) which are strong solutions to 

V / t>Si 

{*,Si,Ui), driven by independent BM's and hence independent. Then by lemma 3.3, the sum 



i=l 

has Laplace transform 

e 

A — ^ JJ e-"'*^-*3W , A > . 
1=1 

2) Arrange the point masses (Tj,!^) of fi in decreasing order cxd > Yi > 12 > . . . of the second 
component. Let (5„ | denote the sequence selected in the proof of lemma 2.4, based on Then 
for every n fixed, we read step 1) conditionally on p,: 

g-A E,6w,Y,>i„ft3" ' I^a) = n e-^^*^-*3W , A>0. 

^ i£lN,Yi>in 
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Under condition ([3]) we have bounds ()22p . and thus as in the proof of lemma 2.4 



oo 



where the hmit is a.s. finite. On the other hand, the collection of points (l^)j>i being PRM on 
(0,00) with intensity [fi^o,{s)ds] ■ iy{dy), ^ always finite by our permanent assumption 

J{y /\ l)i'{dy) < 00. Moreover, s — > ^s,t3i^) is bounded on (ti,t2] for fixed A in virtue of lemma 3.1. 
Thus we have convergence on both sides of (I47p as n — > 00, and arrive at 



(48) 



3) Fix A > and define F : {ti,t2] x (0, 00) [0, 00) by F{s, y) = y ^'s.tglA) . Then we have 
^/'e-/n^,j/)M(rfM?/)'\ = ^(expJ- / \l-e-^^'^yA9{ds,dy)\\ 

^ ^ V I -'(ti.t2]x(o,oo) L J J y 

by [IW 89, p. 42-44]). Taking expectations in (gH]) we have proved ([MD. □ 

3.9 Lemma: a) For all < s < t < 00, there is a law J^,* on [0, 00) with Laplace transform 

(49) ^ e-^"7,,t(dx) = exp(^- J^a{v)^,,t{X)dv^ , A>0 
and the family {Is,t)o<s<t<oo has the same skew convolution property 

(50) It^^tiHt^M * = ^tuta for < ti < t2 < ^3 < 00 . 

with respect to the semigroup (-f^s,t("))o<s<t<oo ^ family {Is,t)o<s<t<oo of lemma 3.6. 
b) The semigroup (e/g^j (•, •))o<s<t<oo of transition probabilities associated to the process 
d^t = -m dt + [ y fxidt, dy) + a{t) dWt , t>0 

J(0,oo) 

is given by 

Js,t{y^ •) = Hs,t{y, ■)* Is,t , < s < t < 00 , y > . 
Proof: By definition of ^i.,t(-) in © and by 3.3 together with (|32p , we have 

^V,t2 ° ^t2,t3 = ^t2,t3 for tl<V<t2<t3, 

thus the skew convolution property (|50p has already been proved in lemma 3.8. This is a) 
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As in lemma 2.2, we prepare independent BM's W^^^ and processes 

i) ^^^^ = F^^\W^^\ s, y) strong solution to equation (*, s, y) of lemma 2.1; by lemma 3.3: 

ii) C^'^^ = F^'^\W^'^\ ^1, s, 0) strong solution to equation ()44p starting in position at time s; combining 
lemmata 2.4 and 3.8 we do have 

c{c?^\cP = o) = I,. 

As strong solutions, both processes ^^^^^^^^ are independent, and their sum ^'•^^ + ^'^^^ provides as 
in lemma 2.2 (this is again the 'branching property') a weak solution to equation (j44|) starting in 
position y at time s. An extension of the Yamada-Watanabe criterion, as in [H 08], gives pathwise 
uniqueness for equation ()44p . Thus assertion b) follows. □ 



3.10 Proof of theorem 1.2: Represent the unique strong solution ^ of equation (o) with initial 
condition y > at time s as sum of two processes ^ = ^^^^ + ^^'^^ driven by independent Brownian 
motions VF(^) and W^^\ where C^^) = F{W'^'^\s,y) is a strong solution to the CIR equation (j24p with 
time dependent coefficients 



(1) 



a{t) - Pit) e« dt + a{t) V(ef V dt^f ^ , t>s, 



and ^(^-^ = F(W^'^\ fi, s,0) a strong solution to (j44p starting in position at time s 



(2) 



m Cj!) dt + [ y f,{dt, dy) + a{t) dwj;^^ , t>s, = 

J{0,oo) 



The branching property as in the last proof (or as in the proof 2.6 of theorem 1.1) combined with 
pathwise uniqueness allows to construct the solution to (o) starting in position y at time s in form 
^ = + and independence of together with lemmata 3.3 and 3.6+3.7 and 3.8+3.9 

identifies the semigroup (Kg^ti-, ■))o<s<t<oo corresponding to the process (o) as 



(51) 



Ks,t{y, 



Hs^tiy,-) * {ls,t*Is,t) , 0<s<t<oo,y>0 



with Laplace transforms 



e Ks,t{y,dx) = exp ( - y^',,t(A) 



dv 



as stated in theorem 1.2. The proof of theorem 1.2 is now complete. 



□ 
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